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. ^ , Abstract 

2^ . We study the single spin asymmetry in the lepton angular distribution of Drell-Yan processes in 

the framework of coUinear factorization. The asymmetry has been studied in the past and different 
results have been obtained. In our study we take an approach different than that used in the existing 
study. We explicitly calculate the transverse-spin dependent part of the differential cross-section with 
suitable parton states. Because the spin is transverse, one has to take multi-parton states for the 

(~| ' purpose. Our result agrees with one of the existing results. A possible reason for the disagreement 

0^. with others is discussed. 

^^ . Single Spin Asyninietry(SSA) can in general exist in high energy scattering with a transversely polar- 

ly-v \ ized hadron. The existence of such an asymmetry implies the existence of helicity-flip interactions and 

o 



nonzero absorptive part of scattering amplitudes. Therefore, the experimental and theoretical study of 
[^ ■ SSA offers a new way to explore the inner-structure of hadrons. Because of its importance, significant 
effort has been devoted to the study of SSA. Reviews about this research field can be found in [1]. 

Theoretical predictions for SSA can be made by using the concept of QCD factorization, if large 
^SJ ■ momentum transfers exist in a process. In the framework of QCD collinear factorization, the nonpertur- 
bative effect of the transversely polarized hadron is factorized into matrix elements of twist-3 operators, 
as pointed out in [21 [3]. In this work we will focus on the SSA in the lepton angular distribution of 
^ ■ Drell-Yan processes. SSA in this case has been studied with collinear factorization in [H [U El [71 El E] , but 
^ . different results have been obtained. The purpose of our work is to solve the discrepancy with a method 
which is different than that employed in the past studies. 
We consider the Drell-Yan process 

hA{PA.s) + hB{PB)^l*{q)+X^r{ki)+i+{k2)+X, (1) 

where Ka is a spin-1/2 hadron with the spin-vector s. hs is unpolarized. We take a light-cone coordinate 
system in which h^ moves in the -|— direction and Hb moves in the — direction. Ha is transversely polarized 
with the spin vector s^ = (0,0, s*i). We employ the Collins-Soper frame to describe the lepton angular 
distribution [To]. In this frame the lepton pair is in rest. We take the spin direction as the direction of 
the X-axis and denote the solid angle of i^ in the frame as Q. We denote the invariant mass of the lepton 
pair as Q and S = {Pa + Pb)"^- We define the following SSA relative to the spin direction: 

^ f da{s±) _ da{-sj_y\ I ( da{sx_) da(-£j_)\ , , 

^ \dQ^dn dQ^dn )l KdQ^dn dQ'^dn J' ^ ^ 



As mentioned, the existing results for this asymmetry are different. Prom [5l[6] the result reads: 

Mm_ sin(2g)sin</. j dxdyTF{x,x)q{y)5{xyS - Q^) 



Q(l + cos2 e) J dxdyq{x)q{y)6{xyS - Q^) 



In the above Tf{x,x) is a twist-3 matrix element of Ha, whose definition will be given later. q{y) is the 
antiquark parton distribution of hs- In |1] the derived A^ has a derivative term of Tp in addition to 
the above expression. Later, the asymmetry has been re-studied in [8l[9]. From [8] Aj\[ is only the half 
of that given in Eq. ([3|) , while the study in |9j confirms the result from [SJ [6] . 

The result for the defined asymmetry has a number of interesting aspects in comparison with the 
asymmetry defined with other differential cross section, which has been studied extensively in |1H [T2t [T3] . 
The asymmetry studied in [11^ \T2[ [13] is at the order of as and has different contributions. In contrast, 
the asymmetry in Eq.(I5]) is at the order of a^ and is predicted in a simple form. Because of these it likely 
provides the best way to access twist-3 matrix elements by measuring Ai\f. Therefore, it is important 
to solve the theoretical discrepancy of this asymmetry. It should be noted that all mentioned results of 
Ai\f are derived with the method of diagram expansion at hadron level, in which one works directly with 
hadron states to evaluate differential cross-sections of hadrons. 

It should be realized that QCD factorizations are general properties of QCD, if they are proven. In 
principle one can derive the factorization for a hadron scattering by replacing hadrons with QCD partons. 
After the replacement one can explicitly calculate the differential cross section of the corresponding parton 
scattering and relevant matrix elements of QCD operators. With the obtained results, one can directly 
derive perturbative coefficient function in the factorization, hence the factorization. It has been started 
in |14t [T5| [TBI [TTl [TS] with this approach to derive the factorization of SSA. In this work we will take this 
approach to study the collinear factorization of Aj\f . 

Before we turn to our calculation of SSA with partonic states, we give the definition of the twist-3 
matrix element appearing in Ai\f. We will use the light-cone coordinate system, in which a vector a'^ is 
expressed as a^" = (a+,a",a_L) = ((a° + a^)/V2,{a'^ - a^)/\/2,a^,a^) and a]_ = {a^f + (a^)^. In the 
system we introduce two light-cone vectors: n^ = (0, 1,0,0) and l^^ = (1,0,0,0). Other notations are: 

gi^^ =gi^'^ -n'^l'^-n^l", e^^;; = e'"^'"'lan,3, e^^^^ = I. (4) 

The definition of the twist-3 matrix element Tp{xi, X2) reads: 

Tf{xi,X2) = -5sS;.|^^e-^^^(^^--'^)^^-*^^-'i^^(P^,sx|V;(yin)7+G+^(y2n)^(0)|P^,sx) (5) 

with s'^ = e^s±,^. The above definition is given in the light-cone gauge n ■ G = 0. In other gauge gauge 
links along the direction n should be supplemented to make the definition gauge invariant. Instead of 
the spin- vector one can also use helicity A to describe the polarization of h^. In this case the general 
forward scattering amplitude like (/iyi(A')|C'|/iyi(A)) with some operator O is a 2 x 2 matrix in the helicity 
space. It is clear that the non-diagonal part corresponds to the forward scattering amplitude of Ha with 
the transverse polarization. 

If we use a transversely polarized single-quark state to replace Ha in the definition of Tp, one will 
always get Tp = because the helicity of a massless quark is conserved in QCD. In order to study SSA 
and its factorization, one has to consider multi-parton states for the replacement [16]. Following [16] we 
consider the state 

\n[X]) = \q{p, A)) + cMpu \)9{P2. \g)[\ = \ + \g\) + • • • , (6) 
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Figure 1: The diagram for Tp{x,x) at one-loop. The black dots denote the insertion of operators used 
to define Tp (xi,X2) in Eq.(l5|). 



with pi+ P2 = P- The state \n[X\) is a superposition of a single quark- and quark-gluon state. It has the 
helicity A. The single quark state must have the same helicity A, while in the quark-gluon state, denoted 
as qg-state, the sum of the quark helicity Xq and the gluon helicity A^ must be A = Ag + A^. The qg-state 
is in the fundamental representation of SU(Nc)-gauge group. We take the momenta as pi = xop and 
P2 = {I — xo)p with pf^ = P^ = (P^, 0,0,0). If we replace /lyi(A) in Eq.(l5|) with the state |n[A]), one 
will find that Tp receives contributions only from the matrix elements of the interference between the 
single quark- and the qg-state, i.e., {q{Xq)\0\q{Xg)g{Xg)) or {q(Xq)g{Xg)\0\q{Xg)) . It is noted that the 
total helicity in the bra- and ket-state is different, but the quark always has the same helicity. In general, 
one can use the state in Eq.([6]) to construct a 2 x 2 spin density matrix (n[A']|0|n[A]) in helicity space 
for a given operator O, as discussed in detail in [16l [171 [18] . By taking corresponding operators, e.g., the 
one used to define Tp, one can obtain from the non-diagonal part the transverse-spin dependent matrix 
element with O. 

With the state \n[X]) instead of Ha in Eq.([5]) one can calculate Tp perturbatively. The tree- level result 
can be found in |161ll7j. At this order the matrix element Tp{x, x) relevant in this work is zero. Tp{x, x) 
becomes nonzero at one- loop. As found in |16l I17j . at one- loop level there is only one diagram giving 
nonzero contribution to Tp{x,x) in the light-cone- or Feynman gauge. The calculation of the diagram 
is straightforward. The contribution has an U.V.- and a collinear divergence. Both are regularized with 
the dimensional regularization as poles of e = 4 — d. After extracting the U.V. pole we have|161ll7j: 



Tp{x,x,fj,) 



gsOis 



Nc{Nl - l)xo^ 



+ ln 



eV 



+ 0{gsa 



sli 



(7) 



where the pole is the collinear divergence with the index c. /i is the renormaliation scale related to the 
U.V. pole, and ^c is that related to the collinear pole. For simplicity we have taken ci = 1 in Eq.(l6]). 
Now we turn to the Drell-Yan process in Eq.(l). The relevant hadronic tensor is defined as: 

W^' = Y. j ^^e^'^^hA{PA,s^)M{PB)\mi''m\X){X\q{x)^^^q{x)\hB{PB),hA{PA,s^)) (8) 



X 



where q is the momentum of the lepton pair. The differential cross section appearing in Eq.Q is related 
to the tensor as: 

(9) 



da 



a 



■ J dUS{q^ - 0") [Kk'i + Kkl^ - h ■ k2gn ^m- 



It should be noted that in the defined distribution in Eq.(9) only the invariant mass Q of the lepton pair is 
fixed and some components of q are integrated, e.g., q±. In the integration one should note that the lepton 



momenta ki^2 depend on the momentum q in the moving frame and on the sohd angle Jl in Collins-Soper 
frame. Because the integration over q with q'^ fixed, the integration can give some soft divergences in the 
small q± region. Therefore one should perform the factorization of the defined differential cross-section in 
Eq.® instead of structure functions of W^'^. Only in the case with other observables which are directly 
related to structure functions, one needs to perform factorizations for these functions. 

We denote the spin-dependent and symmetric part of W^^'^ as W^'^. Only this part will give contri- 
butions to Ajy. For q± ^ the tensor can be decomposed into eight structure functions [121 [20]. These 
structure functions are in general singular with gj_ — )• 0. In principle one can have a part of W^'^ which 
is proportional to 6'^{q±). This part has a simple form: 



W^" 



pu 



pu 



Pa-Q Pb 



+ 5" 



n 



P'b 



Pa-Q PB-q, 



Wq + --- 



(10) 



with Wq (X S'^{q_i). The • • • represent the terms of the eight structure functions which can be found in 
[HIEO]. One expects that these structure functions can have a factorized form. In [HI 121 13] , the trace 
part, i.e., W^^ has been studied with collinear factorization. In [16 1 1171 [T8] . the factorization of the same 
part has been studied with multi-parton states. It should be noted that the factorization of structure 
functions can be different than that of Aj\f. The structure functions are for fixed q±. For A]\f, i.e., for 
the differential cross section in Eq.([9]) qj_ is integrated over. The collinear- and I.R. divergences in Eq.Q 
can appear in a different way than those in structure functions. This will be clearly seen in the results 
obtained with the multi-parton state. 

Having the result for Tf{x,x) of the state \n{p, A)) in Eq.© , we need to calculate the asymmetry 
of the scattering with the same state \n{p, A)) to study the factorization or to check the factorization in 
Eq.Q. For this we consider the scattering: 



n{p,X) + q{p)^^*{q)+X, 



(11) 



with Pa = p and Pb = P- The leading order of SSA here is at ggCis- As we will see, the obtained 
differential cross-section contains collinear divergences. We will show that the collinear divergences can 
be factorized with the collinear divergence in Tf{x,x) of Eq.©. Because Tf(x,x) is at order of gsOig, it 
results in that the perturbative coefficient function is at order of a^ as in Eq. ([3]) . The finite contributions 
in the differential cross-section can be factorized with the tree- level result of Tp(xi, X2) which is at order 
of 5s. Hence, the finite contributions will delivery corrections at order of a^ to the perturbative coefficient 
function. Hence, we only need to find collinear divergences in the differential cross-section at the leading 
but nontrivial order. 
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Figure 2: Classes of diagrams for the contributions to the hadronic tensor with the multi-parton state in 
Eq.®. The black dots are insertion of electromagnetic current operators. The broken line is the cut. 



The contributions to W^^'^ can be classified into three classes of Feynman diagrams for W^'^ , which 
are given in Fig. 2. In each diagram of Fig. 2 there is a cut dividing the diagram into a left- and right part. 
In order to have SSA, a cut should be also exist in the left- or right part. This cut is not drawn in Fig. 2. 
In Class (a) of diagrams, represented with Fig. 2a, there is no parton in the intermediate state. In Class 
(b) of diagrams represented with Fig. 2b, the initial gluon without interactions with partons in the left 
part of the diagram, goes through the cut to interact partons in the right part. The intermediate state 
only contains this gluon. It is clear that contributions from these two classes of diagrams are proportional 
to S'^{q±). Fig. 2c represnts diagrams of Class (c). In these diagrams the intermediate state contains an 
emitted gluon. Hence, in contributions from Class (c) q± can be nonzero. It is possible to have an 
additional class of contributions, which are those diagrams where the initial gluon in the left part in 
Fig. 2c can go without interactions through the cut. But at the leading order, there is no absorptive part 
in the left- or right part. Hence the contributions from this additional class are zero at the order. 

Not all classes of diagrams need to be considered. One can show that the contributions from Class (b) 
are exactly zero. For Fig. 2b we denote the left- and right part as the amplitude Tj^ and T^, respectively. 
These amplitudes are in fact the matrix elements: 

K = {0\J^\q{pi,XMp)), U = {qipAMP)\Jl9{P2,Xa)), (12) 

where we have labeled the helicity of the quark and gluon explicitly. It should be noted that The qq in 
T^ is in color-octet, i.e., the pair carries the same color as that of the gluon in T^. These amplitudes 
can be decomposed into the form factors: 

Ti: = v{p)rTMP, \)Fi{q^). U = (^^ - ^) Hp, \)T''i ■ e{\g)v{p)G^{q^), (13) 

\q-pi q-pj 

where the color index a is the color of the gluon, e(Ag) is the polarization vector of the gluon and q is 
fixed a.s q = pi + p. The form factors Fi and Gi are complex functions of q^ in general. Using these 
expressions one can calculate W^^'^ directly. One easily finds: 

wi;^^{i + \\). (14) 

Since we take the state in Eq.Q as a spin-1/2 system because Ha is with spin-1/2, one always has 
\q\g = —1. Therefore, the contributions from Class (b) are zero. This holds at any order of a^. In fact, 
the conclusion about the contributions from Class (b) is a consequence of the helicity conservation of 
QCD. 

There are many diagrams for contributions from Class (a). Since we are interested in soft divergences, 
we only need to consider those diagrams which contain soft divergences. Those diagrams are given in 
Fig. 3. In each diagram, a short bar cutting a quark propagator means to take the absorptive part of 
the propagator. In fact the short bar implies one or more physical cuts in the amplitude of the left- 
or right part of the diagrams. E.g., when we extend the short bar to the bottom of Fig. 2a, the right 
part represents in fact the scattering 7* — ?• g? ~^ QQi hence a nonzero absorptive part of the amplitude is 
generated. A special care should be taken for the left part of Fig. 3a, because two collinear partons merge 
into a quark propagator. This brings up an ambiguity like 0/0, when two partons are exactly collinear. 
To deal the ambiguity we first take the gluon momentum p2 off-shell by giving it a small — component, 
i.e., P2 = (^2^,^2^,0,0). Then the left part of Fig. 3a becomes: 

rL L) = 9sv{p)rr^^^^rr-X^ • ^^"^(^i) = 9sv{p)i \ . ^y\^^. T^I ■ eT'^u{p,), (15) 

'^•^"•■^^ [Pi + P2Y + i£ 2p2 [pi + P2r + ^£ 




Figure 3: The diagrams of Class(a) contributions. These diagrams contains soft divergences, 
diagrams which do not contain soft divergences are represented with • • •. 
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from the above we can in the last step take p^ = 0. This is equivalent to take the quark propagator as 
the special quark propagator given in |21j . 

The soft divergences in Fig. 3 can easily be worked out. Here we notice that the soft divergences in 
these diagrams are not collinear divergences. They are generated through exchange of a Glauber gluon, 
whose momentum k has the pattern k^^ ~ (A^,A^,A, A) with A — )• 0. The reason for this momentum 
pattern is the physical cuts in amplitudes. With the on-shell conditions from the cuts the +- and — 
component of k must be at order of A^ if k± is at order of A with A — t- 0. For the hadronic tensor we also 
need to add the contributions of conjugated diagrams of Fig.3. We have the divergent part regularized 
with dimensional regularization: 






6a+66 



6e+6/ 



6c+6d 



gsQ^s 



(iV2 _ 1)^ 



4A^2 



■5\p + p- q)V2^s^'' {p - pf^ ( -— ) + 



-QsQ-s 



5 {p + p- q)^/2xoS^'^ {p - p) 



{- 



\iA 



2 
2 



ec 



+ 



gsQ-s 



Nj-l 

Am 



5''{p + p-q)V2xos^''ip-pr +•••' 



m} 



ec 



(16) 



with ec = 4: — d. Its pole represents the Glauber divergence. The terms represented by • • • are finite. We 
have used the notation to denote the symmetric part of a tensor built by two vectors: 



Ait^B"^ = M'B'' + A''B^'. 



(17) 



From Eq. (|16p the sum of the three pieces of contributions are finite. Therefore, we conclude that the 
contributions of Class (a) to the hadronic tensor and hence to the differential cross-section are finite at 
the leading order. 

Now we turn to the contributions from Class (c). At the leading order the contributions come from 
diagrams given in Fig. 4. These contributions to W^^'^ are nonzero and finite for qx_ ^ 0. But when 
we calculate the differential cross-section, a divergence appears after the integration over q^. To find 
the divergence we can in the first step expand those contributions from Fig. 4 m. q± ^ \ with A — t- 
and keep only those contributions which will give the divergence after the integration over q±. From 
Eq.Q the contributions to W^'^ at order of A~^ or higher will give finite contributions to the differential 
cross-section. 
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Figure 4: The diagrams for Class (c) 



By expanding the contribution from each diagram in Fig. 4, we find that all contributions except those 
from Fig.4e are at order of A~^ or higher order. Therefore, we have the divergent part of W^'^ of Class 
(c): 



W^"" 






(N^ - 1)V2^5{1 - y)5{x - xo) 



1 



flU 



1 



^i^'n''^ 



.S^ImsI'} 



[Qi) P-P V -^2 



(gj_)2 ^ 2p-pgj_ 

+ 0(A-i). 



^xop^^s*^} 



(18) 



In the above, the contribution with g^ is at order of A~^ and the contribution from the remaining three 
terms is at order of A"^. Only the terms in the first line of Eq. ljlSp give divergent contributions to 
the differential distribution in Eq.(9) when we integrate over g_L. The terms in the second line do not 
give divergent contributions because of the rotation covariance. The terms at A~^ will only give finite 
contributions. 

Adding all contributions together we have W^^'^ at the leading order: 



W^""" 



QsOis 



Air 
+5{x - xo)6{l -y) 



{Nl - l)^j2Fo\ (n - l)^n''h{l - x)5{l - y)6\q^) 



wq 



/lU 



(gf)2 ^ 2p-pqj^ 



{f^S"} 



^xop'^-s 



,2/2 - (' ■ 'i±p^'''i± + ^^Ip^^s-^ 
{qirp-pK 2 



+ 0(A-i)|+0(g,Q2 



(19) 



wq is a finite constant. Our result of is C/(l)-electromagnetic gauge invariant up to order of A ^. This 
can be checked by noting g^ = p^ + pf^ + qj_ + 0{\'^): 

q^W^'•' = q,W^"' = + O{X-^). (20) 

We notice that the spin- independent part of the hadronic tensor W'^'^ is proportional to 6'^{q±) at the 
leading order of a^, because only the process qq — )• 7* contributes. This is in contrast to the case studied 
in Eq. (jl9p . where one has a contribution with q± ^ 0. 

Substituting the result in Eq. (|19|) into Eq.Q and performing the integration over g_|_, we find 

2_\1 SjxQS - Q2 ^ 
ej\ 1287r2g3 



da{s±) 
dQ^dn 



-{N^ - 1)2:0 ^/2i;; 



■ sin(20) sine/) + • • • 



(21) 



where the divergence is regularized with ec = 4 — d and • • • denote finite contributions. The finite 
contributions are from the term with wq and those at order 0{\~^) in Eq. (jl9p . We notice here that the 
upper limit of q± in the integration is finite by the energy- momentum conservation. In Eq. (j2ip we take 
the spin-direction in the x-direction, and the moving direction of i~ in the Collins-Soper frame is given 
by (sin 9 cos cp, sin sin (p, cos 9) . 

The divergence comes from the g^-region with g^ ~ A — t- 0. In this region one easily finds that 
the gluon exchanged at the bottom of Fig.4e is collinear to the incoming gluon. Its momentum k scales 
like k^^ ~ 0(1, A^, A, A) with A — )• 0. Therefore, the divergence is a collinear one. With the result 
of Tf{x,x) given in Eq.d?]) and the leading result for an antiquark distribution in an antiquark, i.e., 
q{y) = (5(1 — y) + 0{as), we can derive the factorized result for the differential cross-section: 

da{s±) o? 

and the asymmetry: 



sm{29)sm(l) ( dxdyTF{x,x)q{y)5{xyS - Q^), (22) 



sin(2e)sin</) j dxdyTF{x,x)q{y)5{xyS - Q^) 

^^ ~ ~o7vi~~^ 2^ 7 ■ (^^^ 

2(^(1 + cos 9) / dxdyq{x)q{y)6{xyS - Q^) 

This is our main result. In our results we have taken the electric charge fraction Cq of quark as 1. One 
can easily generalize the above results to any flavor. The finite contribution in Eq. (j2ip will be factorize 
with Tf{xi,X2) at tree-level and give a correction at order of Ug to the above A^. Here, we have a case 
that from an partonic observable at a given order of Og the extracted perturbative coefficient functions 
can be at different orders. This is the nontrivial order-mixing as observed in [17J. 

Before we make a comparison with existing results, we show in the below that the differential cross- 
section is indeed factorized with Tp(x,x). For this we write down explicitly the contribution of Fig.4e: 



M 



^' rr''T^K{Pi-q?)5{v--q--k-' 



4e 



2iVc(27r)4-' ^^^' ^J > Kf -i ^k+{p-k)^{p2-ky 



Tr 



T''u{pi)u{p)T jp'j ■ (p-k) 



l^v{p)v{p)^p^T''-i ■ (pi - g)7'" I 



{{P2 + kY'e" + {-2k + p2) ■ egPP' + (-2p2 + kye"') , (24) 

where k is the momentum of the gluon crossing the cut with k'^ = p'^ — q^ and k± = —q±. We have 
performed the integration over k~ with the on-shell condition of the gluon. This gives k~ = k'^/{2k~^). 
For finding the collinear divergence one can safely neglect k" in the (5-function for — components of 
momenta. The product in [• • •] in the above is a matrix with color- and Dirac indices. This matrix can 
be expanded with 7'^and 7^75 for Dirac indices and with the color matrix T" for color indices. We can 
write the product as: 



/ 



abc 



T^u{pi)u{p)T\-f {p - k) 



jabdr^ (^J,drj.arj.h\^ rpC (^^- A\ + B p) 



A+ = u{p)-fp-f ■{p-k)j+u{pi). (25) 



In the above we have explicitly written out the structure 7 for Dirac indices, the remaining terms have 
been denoted as Bp. If we only keep the term with A'^ and take 7pj in Eq. (p^ as 7 • lup-^, i.e.. 



M 



= - 2N,{2'Kf '^^^^^^ ~ ^^'^'^^^" " "^'^^ {T'l-l''v{p)v{p)rT^l ■ {Vi - qhn 



i{p2 + k)P^eP + i-2k + p2) ■ egPP^ + {-2p2 + k)PeP^)\ + • • • , (26) 

where • • • represent the remaining contributions from Bp in Eq. ()25p or 7p^ in Eq. (j24p with pi =_L. We 
have: 



^^' = -^-^{Nl-l)V2^o5{l-y)5{x-xo) 



^ ~ i^-" ^ —}^—^^r^{i^Z''} 



{qj_Y 2p-pqj^ 

+0{X-^)+0{X-^). (27) 



The contributions in the first hue come only from the term with A'^ in Eq. (j25p and with 7^^ in Eq.(p4 
as 7 • lup-^ . These contributions are exactly those in the first line of Eq. (jlSp , which generate the divergent 
contributions in the differential cross-section given in Eq. (j2ip . The terms represented by 0{X^^) do 
not give divergent contributions with the similar reason as discussed after Eq. (ll8p . Now we note that 
the integral over k± with the part in {• • •} in Eq. (j26p multiplied with a factor kj_ is proportional the 
contribution to Tf{x,x) given in Fig.l. The factor k^ can come from the first line of Eq.(|26p. or from 
the leptonic tensor in Eq.Q, when we express the lepton momenta with the lepton momentum in the 
Collins-Soper frame and the momentum q. It should be noticed that k± in Eq. ()26p is just —q±. Therefore, 
the obtained divergence in the spin-dependent part of the differential cross-section is exactly factorized 
with Tp{x,x). With the above discussion, it is also clear that for the differential cross-section in Eq.Q 
one should perform the factorization only after the integration over g^. 

Our result disagrees with those in [H O [H [9]. Aj^ obtained here is only half of those derived in 
[U El [9]. Although we have used a method different than that in the existing studies, our result agrees 
with that given in [8], where the contributions factorized with chirality-odd operators are also given. In 
[9], the issue of gauge invariance in the problem has been emphasized. In our study, gauge invariance is 
explicitly kept because we take on-shell parton states to calculate the hadronic tensor. Our result is also 
C/(l)-gauge invariant as checked in Eq.(20). 

We can explore the reason for the discrepancy. In [5l EJ [9] one factorizes the hadronic tensor in the 
first step. Then in the next step one uses the factorized tensor to calculate the differential cross-section 
by performing the integration over g_L indicated by Eq.([9]), where the factorized tensor seems proportional 
to Si^{q±)- This procedure is only correct if the integration over q_\_ does not generate soft divergences. 
However, the integration does generate a collinear divergence as we have seen here. This is the main 
reason for the discrepancy. To see the discrepancy more clearly, we re-arrange our result in Eq. (|19p as: 



^p. = -'-^{Nl-l)V2^o5{l-y)5{x-x^] 
1 



■ nn'l V / 



fiu s ■ q± _| u} 
[qiY V P-P 



2p ■ pqj_ 
+ •••, (28) 



where • • • stand for terms which are irrelevant here or at higher order of a^- In [5l [6l [9] the factorized 
tensor used to calculate An has only the tensor structure which is equivalent to the part given in the 
first line in the above. If we only use this part to calculate A^, we obtain the same factorized result. It 
is noted here that the soft divergence of this part is proportional to S{q'^) ~ S'^{q±) if we take this part 
as a distribution of q'j_. That is, this part may be factorized before the integration over qj_. But, the part 
in the second line also give a divergent contribution to ^^v and it results in the discrepancy. 



To summarize: We have studied SSA in the lepton angular distribution of Drell-Yan processes. The 
asymmetry has been studied with the diagram expansion at hadron level. Different results have been 
obtained. In this work, we take a different approach. We calculate the transverse-spin dependent part 
of the differential cross-section with suitable parton states. Because the spin is transverse, one has to 
take multi-parton states for the purpose. Our result agrees with that in [8], but disagrees with those in 
[H El El [9]. A possible reason for this has been discussed. It should be emphasized that the studied SSA 
is very interesting, because it is at order of a^ and its prediction takes a simple form. Measuring such a 
SSA likely provides the best way to access the twist-3 matrix element. Finally, we notice that there is 
also a corresponding SSA starting at order of q^ in semi-inclusive DIS and its prediction takes a simple 
form. The work about this will appear elsewhere. 
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